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This  paper  presents  a  general  model  for  air-side  feed  stream  contamination  that  has  the  capability  of 
simulating  both  transient  and  steady-state  performance  of  a  PEM  fuel  cell  in  the  presence  of  air-side 
feed  stream  impurities.  The  model  is  developed  based  on  the  oxygen  reduction  reaction  mechanism, 
contaminant  surface  adsorption/desorption,  and  electrochemical  reaction  kinetics.  The  model  is  then 
applied  to  the  study  of  air-side  toluene  contamination.  Experimental  data  for  toluene  contamination  at 
four  current  densities  (0.2, 0.5, 0.75  and  1.0  A  cm-2 )  and  three  contamination  levels  ( 1, 5  and  10  ppm)  were 
used  to  validate  the  model.  In  addition,  it  is  expected  that,  with  parameter  adjustment,  this  model  can 
also  be  used  to  predict  performance  degradation  caused  by  other  air  impurities  such  as  nitrogen  oxides 
(NOx)  and  sulfur  oxides  (SOx). 
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1.  Introduction 

Proton  exchange  membrane  fuel  cell  (PEMFC)  contamination 
caused  by  impurities  is  an  important  issue  in  fuel  cell  operation  and 
applications.  Contaminants  can  be  classified  into  two  categories 
according  to  the  source:  external  and  hardware.  External  impuri¬ 
ties  enter  the  fuel  cell  mainly  through  the  fuel  or  air  feed  stream.  The 
major  impurities  in  hydrogen  fuel  produced  from  reformate  include 
CO,  C02,  H2S,  NH3,  organic  sulfur-carbon,  and  carbon-hydrogen 
compounds,  and  those  in  air  include  pollutants  such  as  NO*,  SO*, 
NH3,  03  and  small  organic  compounds.  The  impurities  introduced 
through  hardware  sources  contain  metal  ions  such  as  Fe3+,  Cu2+, 
Ni2+,  and  Cr3+  from  gas  diffusion  layers,  bipolar  plates,  and  system 
components,  Na,  Ca,  and  S  from  membranes,  and  Si  from  gaskets 
and  cooling  fluid. 

It  has  been  demonstrated  that  even  trace  amounts  of  impurities 
in  the  fuel  or  air  stream  or  the  fuel  cell  system  components  can 
severely  poison  the  anode,  membrane,  and  cathode,  particularly  at 
low-temperature  operation.  The  contaminants  can:  strongly  or  irre¬ 
versibly  adsorb  on  the  catalyst  surface,  blocking  the  reaction  sites; 
enter  the  membrane,  reducing  proton  conductivity;  and  change 
the  hydrophobicity/hydrophilicity  of  the  catalyst  layer  interface, 
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affecting  mass  transportation.  As  a  result,  contaminants  cause  per¬ 
formance  degradation  and  sometimes  permanent  damage  to  the 
membrane  electrode  assemblies  (MEAs)  [1-3].  To  address  this 
issue,  a  great  deal  of  effort  has  been  put  into  identifying  the  poten¬ 
tial  impacts  of  contamination,  understanding  the  contamination 
mechanisms,  and  developing  mitigation  strategies. 

At  present,  contamination  research  concentrates  on  three  major 
areas:  (1)  experimental  impact  observation  and  validation;  (2) 
theoretical  modeling  of  contamination  to  provide  a  fundamen¬ 
tal  understanding  of  the  mechanisms  and  to  develop  predictive 
tools;  and  (3)  contamination  mitigation  strategy  development. 
With  respect  to  contamination  impact  observation,  the  effects  on 
fuel  cell  performance  of  impurities  such  as  CO,  C02,  H2S,  and  NH3 
in  the  fuel  feed  [4-20]  and  SOx,  NOx,  H2S,  and  NH3  [16,21-24]  in 
the  air  stream  have  been  extensively  examined.  A  study  of  battle¬ 
field  air  impurities  such  as  benzene,  propane,  HCN,  CNCL,  sarin,  and 
sulfur  mustard  has  also  been  reported  [25  ].  In  terms  of  model  devel¬ 
opment  for  predicting  fuel  cell  contamination,  numerous  studies 
on  fuel  cell  contamination  have  been  conducted,  especially  for  CO 
contamination  [26-37].  For  fuel-side  H2S  contamination,  a  kinetic 
model  has  been  developed  to  study  both  transient  and  steady  fuel 
cell  performance  [38].  But  to  the  best  of  our  knowledge  a  compre¬ 
hensive  predictive  model  for  air-side  feed  stream  contamination  is 
not  yet  available  in  the  literature,  other  than  a  preliminary  adsorp¬ 
tion  model  for  N02  [39].  Since  air  pollutants  have  a  large  effect  on 
fuel  cell  applications,  for  example,  portable  power  units  for  indoor 
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Nomenclature 

Co2  oxygen  concentration  in  the  catalyst  layer 

CH+  proton  concentration  in  the  CCL 

Cp  concentration  of  species  P 

CP/  concentration  of  species  P 

Cp  contaminant  concentration  in  the  cathode  channel 
^o2,air  bulk  oxygen  diffusivity 
Dg  air  effective  oxygen  diffusivity 
F  Faraday’s  constant 

Jcen  cell  current  density 

ka3f  and  ka3h  forward  and  backward  reaction  rates  of  the 
reaction  (a-iii) 

1<ap3f  and  kap3b  forward  and  backward  reaction  rates  of  the 
reaction  (a-vi) 

kjf  and  /<jb  forward  and  backward  reaction  rates  of  the  jth 
reaction  (j  =  1,  2, . . .,  5) 

kpjf  and  /<pjb  forward  and  backward  reaction  rates  of  the  reac¬ 
tions  (a-iv),  (a-v),  (a-vii)  and  (a-viii)  (j  =  1,  2, 4,  5) 
f'CGDL  thickness  of  the  cathode  gas  diffusion  layer 
n  number  of  the  Pt  sites  occupied  by  contaminant  P 

rij  electron  transfer  number  for  the  corresponding 

individual  electrochemical  half-reaction  (j  =  1, 2, . . ., 
5) 

na3  electron  transfer  number 

nap 3  electron  transfer  number  of  the  reaction  (a-vi) 

no2  electron  transfer  number  of  the  ORR 

Po2  oxygen  partial  pressure 

Pcgc  gas  pressure  in  gas  channel 

Pvapor  saturated  water  vapor  pressure 

q  electron  transfer  number  of  the  reaction  (a-vii) 

R  ideal  gas  constant 

Ro  fuel  cell  internal  resistance  (membrane  resistance) 

t  time 

T  temperature 

Vceii  cell  voltage 

V°  open  circuit  voltage  (OCV) 

z  electron  transfer  number  of  the  reaction  (a-viii) 

Greek  symbols 

otj  electron  transfer  coefficient  for  the  corresponding 

individual  electrochemical  half-reaction,  j  =  1, 2, . . 

5 

aa3  electron  transfer  coefficient  of  the  electrochemical 

half-reaction 

otap3y  ap4  and  ap5  electron  transfer  coefficients  of  reactions 
(a-vi),  (a-vii)  and  (a-viii) 
r  total  surface 

yc  ratio  of  the  active  surface  to  the  geometric  surface 

of  the  CCL 

£cgdl  porosity  of  the  cathode  gas  diffusion  layer 
r] a  anode  overpotential 

q c  cathode  overpotential 

0pt,  #pt-o2»  #pt— o2h>  #Ptn-P  and  0ptn-p'  surface  coverages  of 
corresponding  species 

0pt ,  #pt_o2  and  0pt_o2H  surface  coverages  of  corresponding 
species  without  contaminant 


and  underground  applications  and  fleet  vehicles  for  outdoor  appli¬ 
cations,  developing  a  predictive  model  for  air-side  contamination 
is  necessary  in  order  to  evaluate  fuel  cell  performance  under  the 
above-mentioned  conditions. 

In  this  paper,  we  present  our  work  on  the  development  of  a  gen¬ 
eral  model  for  fuel  cell  air-side  feed  stream  contamination.  This 
is  a  kinetic  model  that  is  mainly  based  on  the  proposed  oxygen 
reduction  mechanisms,  the  contaminant’s  surface  adsorption,  and 
the  electrode  reaction  mechanism.  The  model  has  the  capability  to 
simulate  both  transient  and  steady-state  cell  performance  in  the 
presence  of  an  air  contaminant.  The  model  is  further  validated  by 
our  experimental  data  collected  in  the  presence  of  toluene  in  a  fuel 
cell  air-side  feed  stream.  These  experimental  data  were  obtained 
at  different  current  densities  and  several  toluene  contamination 
levels.  The  model  is  further  employed  to  examine  the  toluene  con¬ 
tamination  effect  on  cell  performance  at  low  concentration  levels. 

2.  Model  development 

2.1.  Oxygen  reduction  reaction 

The  oxygen  reduction  reaction  (ORR)  is  a  multi-electron  multi- 
step  reaction.  For  the  ORR  on  a  platinum  (Pt)  surface,  two  Tafel 
regions  are  generally  observed  in  both  acidic  and  alkaline  solu¬ 
tions.  At  low  current  densities  (high  electrode  potential),  a  Tafel 
slope  of  -60  mV  dec-1  can  be  observed,  whereas  high  current  den¬ 
sities  (low  electrode  potential)  yield  a  Tafel  slope  of  -120  mV  dec-1 
[40].  The  difference  in  Tafel  slopes  is  attributable  to  partial  cov¬ 
erage  of  the  Pt  surface  by  the  oxygen-containing  intermediates  O 
and  OH.  Experimental  study  indicates  that  the  reaction  rate  is  first 
order  with  respect  to  oxygen.  Extensive  research  of  the  ORR  on  Pt 
surfaces  has  been  conducted  in  the  past  40  or  so  years  and  a  con¬ 
siderable  number  of  ORR  mechanisms  have  been  proposed  [41  -44]. 
Yet  despite  numerous  studies  of  the  ORR,  the  detailed  mechanism 
remains  elusive.  Damjanovic  and  co-workers  [40,45,46]  proposed 
that  the  charge  transfer  to  the  adsorbed  oxygen  molecule,  with  or 
without  simultaneous  proton  transfer,  is  the  rate-determining  step 
(associative  model): 

M  +  02^M-02  (1) 

M-02  +  H+  +  e-  -*  M-02H(rds)  (2) 

M-02H  +  3H+  +  3e-  ->  2H20  +  M  (3) 

where  M  represents  a  metal  such  as  Pt.  Yeager  et  al.  [47]  proposed 
that  the  first  step  involves  dissociative  chemisorption  of  the  02 
molecule,  which  occurs  simultaneously  with  the  charge  transfer 
(dissociative  model): 

Pt+  1o2  +e-  ->•  Pt-cr  (4) 

In  the  search  to  understand  ORR  mechanisms,  several  model- 
based  impedance  spectroscopy  studies  were  reported  that 
employed  either  an  associative  kinetic  model  or  a  dissociative  one 
[48-50].  In  another  kinetic  model  study  [51],  both  associative  and 
dissociative  mechanisms  were  analyzed. 

In  recent  years,  many  theoretical  explorations  of  the  ORR  mech¬ 
anism  using  quantum  mechanic  methods  have  been  reported, 
as  shown  in  the  review  articles  [52,53].  These  studies  supplied 
information  on  each  elementary  step,  such  as  activation  ener¬ 
gies,  reaction  energies,  and  reversible  potentials.  Anderson  and 
co-workers  [54,55]  investigated  the  activation  barrier  for  each  of 
the  following  electron  transfer  steps: 

Pt-02  +  H+(aq)  +  e"  (U)  -*  Pt-OOH 

Pt-OOH  +  H+(aq)  +  e"  (U)  -*  Pt-(OHOH) 


(5) 

(6) 
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Pt-OHOH  +  H+(aq)  +  e”  (U)  -*  Pt-OH  +  H20  (7) 

Pt-OH  +  H+(aq)  +  e”  (U)  ->  Pt-OH2  (8) 

Their  results  showed  that  the  first  electron  transfer  has  a  smaller 
activation  energy  barrier  than  that  of  oxygen  dissociation,  and  that 
OOH  easily  dissociates  once  formed  after  the  first  electron  transfer. 

Norskov  et  al.  [56,57]  examined  the  thermodynamic  properties 
of  both  associative  and  dissociative  mechanisms: 


\  02+*  -►  o*  (9) 

0*  +  H+  +  e"^  HO*  (10) 

HO*  +  H+  +  e-^  H20+*  (11) 

where  *  denotes  a  site  on  the  surface. 

02  +  *^02*  (12) 

02*  +  (H+  +  e_)  -*  H02*  (13) 

H02*  +  (H+  +  e_)  ->  H20  +  O*  (14) 

0*  +  (H+  +  e-)  -*  HO*  (15) 

HO*  +  (H+  +  e-)  ->  H20  +  *  (16) 


The  results  illustrated  that  both  mechanisms  are  possible  and  are 
controlled  by  electrode  potential.  At  high  electrode  potential,  the 
activation  energy  barrier  for  02  dissociation  will  increase;  thus,  the 
associative  mechanism  is  dominant.  Jacob’s  study  [58]  also  demon¬ 
strated  that  both  the  associative  and  the  dissociative  mechanism 
may  occur  simultaneously,  with  a  higher  ratio  of  occurrence  for 
the  associative  pathway.  In  another  theoretical  study  [59],  Walch 
et  al.  formulated  a  suitable  mechanism  table  for  the  prominent  ORR 
pathway.  This  mechanism  includes  chemical,  electrochemical,  and 


side  reactions  as  follows: 

Chemical  reaction: 

02(ads)  **  O(ads)  +  0(ac[s)  (17) 

°2H(ads)^  °(ads)  +  °H(ads)  O8) 

H20(ads)^  H(ads)  +  HO(ads)  (19) 

OH(ads)  +  OH(ads)  O(ads)  +  ^20(ads)  (20) 

OH(ads)  +  0(ads)  0(ads)  +  OH(ads)  (21) 

Electrochemical  reaction: 

°2(ads)  +  H+(aq)  +  e_  **  °2H(ads)  (22) 

0H(ads)  +  H+(aq)  +  e_^  H20(ads)  (23) 

°(ads)  +  H+(aq)  +  e_  0H(ads)  (24) 

H20(ads)^  0(ads)  +  2H+(aq)  +  2e_  (25) 

Side  reaction: 

Pt  Pt2+(aq)  +  2e-  (26) 

Pt2+(aq)  +  H20  Pto  +  2H+(aq)  (27) 

Pt2+(aq)  +  H2^Pt  +  2H+(aq)  (28) 


2.2.  Proposed  general  cathode  contamination  model 

In  our  attempt  to  develop  a  cathode  contamination  model,  we 
start  with  a  simple  ORR  mechanism.  Since  the  associative  model 


has  substantial  experimental  and  theoretical  support,  we  adopt  it 
and  propose  the  following  mechanism: 


fcif 

Pt  +  02^Pt-02 

fcib 

k2fexp(-(a2n2Fr)c/RT)) 

Pt-02  +  H+  +  e-  ^  Pt-02H 

k2  b  exp(((l  -a2  )n2Fr]c)/ RT) 
k3fexp(-{a3n3Frjc/RT)) 

Pt-02H  +  e-  +  H+  H20  +  Pt-0 

k3b  exp(((  1  -a3  )n3Fr]c )/ RT) 
k4f  exp( -(a4n4  rjcF/RT )) 

Pt-0  +  H+  +  e-  Pt-OH 

k4h  exp(((  1  -a4 )n4 r]cF)/ RT) 
k5fexp(-(a5n5r]cF/RT )) 

Pt-OH +  H+  +  e“  <=►  Pt  +  H20 

k5h  exp(((  1  -a5  )n5  rjcF)/RT ) 


(i) 


(ii) 


(hi) 


(iv) 


(V) 


where  Pt  is  the  active  site  for  oxygen  adsorption,  which  may  contain 
more  than  one  Pt  site,  /<jf  and  /<jb  are  the  forward  and  backward  reac¬ 
tion  rates  of  the  jth  reaction  (j=  1, 2, . . .,  5),  otj  is  the  electron  transfer 
coefficient  for  the  corresponding  individual  electrochemical  half¬ 
reaction,  nj  is  the  electron  transfer  number  for  the  corresponding 
individual  electrochemical  half-reaction,  pc  is  the  cathode  over¬ 
potential,  and  F,  R  and  T  are  Faraday’s  constant,  the  universal  gas 
constant,  and  cell  temperature,  respectively. 

In  our  model,  we  further  propose  that  the  first  electron  trans¬ 
fer  (reaction  (ii))  is  the  rate-determining  step.  In  this  case,  we  can 
further  simplify  the  model  to 


kif 

Pt  +  02^Pt-02  (a-i) 

fclb 


k2f  exp( ~(a2  n2Fr]c/RT)) 

Pt-02  +  H+  +  e_  <=►  Pt-02H  (a-ii) 

k2b  exp(((l  -a2  )n2Frjc )/ RT) 


ka3{exp{-(ua3na3Fr]c/RT)) 

Pt-02H  +  3e-  +  3H+  <=t  2H20  +  Pt  (a-iii) 

ka3b  exP(((l-aa3  )na3Fr]c)/RT) 

where  /<a3f  and  ka3b  are  the  forward  and  backward  reaction  rates 
of  the  reaction  (a-iii),  and  aa3  and  na3  are  the  electron  transfer 
coefficient  and  electron  transfer  number  of  the  electrochemical 
half-reaction,  respectively. 

Reaction  (a-ii)  is  the  rate-determining  step.  The  advantage  of 
this  simplified  model  is  that  we  can  use  one  intermediate  cover¬ 
age  Pt-OOH  as  the  total  oxygen-containing  intermediates  coverage, 
which  also  includes  Pt-0  and  Pt-OH. 

For  a  general  air-side  contaminant  P,  we  consider  the  following 
possible  electrode  and  surface  reactions: 

fcpif 

nPt  +  P  ?=t  Ptn-P  (a-iv) 

kplb 


kp2f 

P  +  nPt-02  Ptn-P  +  n02  (a-v) 

kp2b 


kap3fexP(-(aap3nap3FTlc/RT)) 

P  +  nPt-02H  +  3ne- +3nH+  <=►  2nH20  +  Pt„-P  (a-vi) 

kap3b  exP(((l -“ap3  )nap3Frlc)/RT) 


kp4f  exp(  ( ap4  qFr]c/RT)) 

Ptn-P  +  mH20  +2  Ptn-F  +  qH+ +  qe_  (a-vii) 

kp4h  exp(-((l -ap4)qFrjc)/RT) 


kp5fexp(ap5zFrjc/RT) 

Ptn-P  + 1 H20  nPt  +  P"  +  zH+  +  ze-  (a-viii) 

kp5b  exp(-((l -«p5 )zFrjc)/RT) 


where  Ptn  is  the  adsorption  site  for  P. 
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Based  on  the  proposed  ORR  and  contaminant  mechanisms, 
the  Pt  surface  is  covered  by  the  following  five  species:  Pt,  Pt-02, 
Pt-02H,  Ptn-P  and  Ptn-P,  where  P  is  an  oxidized  form  of  P.  Thus 
for  the  cathode  catalyst  surface  coverage,  we  have 


1  =  0Pt  +  $Pt-02  +  #Pt-02H  +  #Ptn-P  +  0ptn-P'  (29) 

where  0Pt,  0Pt_o2,  #pt-o2H>  0ptn-p  and  0ptn-p'  are  the  surface  coverage 
of  the  above  five  species,  respectively. 

If  we  assume  reaction  (a-i)  is  fast  and  always  at  its  equilibrium, 
we  have 


/<ifCo26>ptr  =  k,b0Pt-o2r  (30) 

where  V  is  the  total  surface  and  Co2  is  the  oxygen  concentration  in 
the  cathode  catalyst  layer  (CCL).  We  will  have 


$pt-o2  = 


kifO)2 

^lb 


Opt 


(31) 


Assuming  reaction  (a-iii)  is  always  at  its  equilibrium, 


■  (32) 


where  CH+  is  the  proton  concentration  in  the  CCL  Solving  the  above 
equation  for  0Pt-o2H  results  in 


$pt— o2h  = 


fcq3b 

^a3f^+ 


exp 


( 


na3r]cF 

RT 


(33) 


Substituting  Eqs.  (31 )  and  (33 )  into  Eq.  (29)  and  solving  for  0Pt  yields 


where  0{Jt,  0pt_O2*  and  0pt_O2H  are  the  surface  coverages  of  the  cor¬ 
responding  species  without  contaminant. 

Based  on  the  proposed  reaction  mechanism,  the  changes  in  the 
surface  coverages  of  Ptn-P  and  Ptn-P  are  given  as 

- =  kplfCp0pt  -  fcplb#Ptn-P  +  /<p2fCp0pt— q2  -  /<p2b^Ptn-PC02 


+  /<ap3fCp0^t_o2HCH+  eXP  - 


0^ap3^ap3rlcF  \ 

-RT  J 


u  pvn  f  (1  aap3)nap3*lcF  ,  „ 

'(ap3b  exP  [  rj  ]  Optn-P 


-kp  4fexP 


ap4qricF 

RT 


&  Ptn-P 


+Vtc!,.xp  (dlzjgw 


RT 


Ptn-P' 


-kpspexp  (  )  6»Pt„_p 


+  kp5bCFC^+  exp 


-(1  -  aps)zr]cF 


RT 


^plb  +  ^p2b^o2  +  ^ap3b  ^XP 


+kp4fexp  (ap4^cF  )  +  kp5f exp  f°^T\ 


(1  —  OCap3)^ap3P!cF  \ 

rt  y 

^Ptn-P 


0Pt  —  A(1  —  0ptn-p  —  0pt(1_pO 
where 


A  = 


1 


1  + 


fclfc02 

fclb 


+ 


^a3b 

fca3f^+ 


exP  ( 


na3hc  F\ 
RT  ) 


Let 


and 


(34) 


(35) 


(36) 


kptfCpA  +  kp2fCpB  +  kap3fCpCH+C  exp 
-)  +kP5bCp-C^exp 

-(l-Otp5)zi/cF^]  „  „  y. 

- ^ -  I  (1  -  ^Ptn-P  -  ^Ptn-P') 

+  VbCH-  eXP  (  — - ^Q’kF  )  Op 


JPtn~P' 


(43) 


c  = 


Fg3b 

^a3fQj 


exp 


H+ 


fna3r]cF 
V  RT 


(37) 


Then  the  surface  coverages  of  Pt-02  and  Pt-02H  are  given  as 
follows: 


$Pt-02  =B(  1  -  0ptn-p  -  0ptn-pO 

(38) 

and 

#Pt-02H  =  C(1  -  0ptn-p  -  0ptn-p') 

(39) 

If  there  is  no  contaminant  present  in  the  cathode  catalyst  layer 
(CCL),  the  surface  coverages  of  Pt,  Pt-02  and  Pt-02H  are  only  deter¬ 
mined  by  the  oxygen  reduction  reactions  ((a-i),  (a-ii)  and  (a-iii))  and 
are  given  as 

q0  _  1 

Pt  1  +  (kifCo2/kib)  +  (ka3b/ka3fC^+)exp(na3r]cF/RT) 

(40) 

q0  k\f  r  ao 

ypt-o2  “  klbC°2t'pt 

(41) 

nO  _  ka3b  exD(n^0cF\nO 

Pt-°2"~  ka3fCiH+  expl  RT  )°n 

(42) 

and 

(44) 

with  the  following  initial  conditions 

0Ptn-p(t  =  0)  =  0  (45) 

and 


0ptn-At  =  0)  =  0  (46) 

where  t  denotes  time  and  Cp  and  Cp-  are  the  concentration  of  species 
P  and  P',  respectively.  The  cathode  current  density  controlled  by 
forward  and  backward  rates  of  reaction  (a-ii)  [60]  is  calculated  by 


Ic  =  n0  2Fyc 


k2f0pt-o2CH+  exp 


(  «2n2'?cF\ 

{  RT  J 


-  k2b0pt-o2H  exP 


/(l  -  ot2)n2r]cF 
\  RT 


(47) 


where  no2  is  the  electron  transfer  number  of  the  ORR  and  yc  is  the 
ratio  of  the  active  surface  to  the  geometric  surface  of  the  CCL. 
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When  there  is  no  contamination,  no  current,  and  overpotential 
(Ic  =  Pc  =  0),  /<2b  can  be  expressed  from  Eq.  (47)  as 


k2b  =  fe2ffl0Pt  ~  CH+ 

^Pt— 02  H 

Substituting  Eqs.  (36)-(42)  and  (48)  into  Eq.  (47)  results  in 


Ic  =  no2Fyck2fl^Co2ACH+  exp  (- 


ot2n2rjcF\ 


1  -  exp 


RT  J 

(1  -  $Ptn-P  -  ^Ptn-P') 


The  fuel  cell  voltage  is  estimated  by 
Vcell  =  -  rja  -  T]c  -  R0IC 


(48) 


(49) 


(50) 


where  V°  is  the  open  circuit  voltage  (OCV)  of  the  fuel  cell,  is  the 
anode  overpotential,  and  R0  is  the  fuel  cell  internal  resistance.  The 
OCV  is  given  by  [61] 


V° 


4.1868 


‘70650  +  8Tlog(T)  -  92.841’ 
.  2F  . 


(51) 


Anode  over-potential  is  generally  small  compared  to  cathode  over¬ 
potential  and  can  be  neglected: 


r)a  «  0 


(52) 


Toluene  adsorption  configuration  is  a  complicated  issue, 
strongly  dependent  on  toluene  concentration,  temperature,  and 
electrode  potential.  Vertical  (a)  and  flat  (tt)  configurations  are  pos¬ 
sible.  Initially,  we  will  consider  only  toluene  adsorption  as  the 
dominant  contamination  mechanism.  We  start  with  the  simplest 
single-site  adsorption  model. 

In  our  model,  oxygen  pressure  in  the  cathode  catalyst  layer, 
which  is  different  from  the  pressure  in  the  channel  due  to  the  pres¬ 
sure  drop  across  the  gas  diffusion  layer,  is  estimated  by  Fields  law 
as 


Po2  =  0.21  (Pccc  -  Pvapor)  -  RT  (54) 

4FD02,air 

where  p q2,  pCGC  and  plfpor  are  the  oxygen  partial  pressure,  gas 
pressure  in  the  gas  channel,  and  saturated  water  vapor  pressure, 
respectively,  Lcgdl  is  the  thickness  of  the  cathode  gas  diffusion  layer, 
and  Df  air  is  the  effective  oxygen  diffusivity  in  air. 

Pvapor  =  exp  (ll  .6832  -  ^g^)  (55) 

The  effective  oxygen  diffusivity  is  estimated  by  the  Bruggeman 
relation, 

D02,air  =  (fiCGDL)1-5D02,air  (56) 


Therefore,  Eq.  (50)  can  be  rewritten  as 
VceU  =  V°  -ric-Rok  (53) 

2.3.  Toluene  contamination  model 

Several  studies  of  toluene  adsorption  on  Pt  have  been  reported 
[62-67].  These  demonstrate  (1)  that  in  the  gas  phase,  toluene  and 
benzene  adsorb  first  as  a  TT-complex  (benzene  ring  parallel  to  the 
surface)  with  a  weaker  electron  transfer  to  Pt  as  compared  with 
ethane  and  propane  [64] ;  (2)  that  an  excess  of  the  toluene  can  trans¬ 
form  the  TT-complex  into  a  cr-complex;  and  (3)  that  the  substituent 
groups  play  an  important  role  in  the  ordering  of  the  overlayers  but 
have  less  effect  on  the  adsorption  strength  [62].  While  adsorbed 
aromatic  compounds  are  in  an  electrochemical  environment,  their 
molecular  orientation  and  packing  density  depend  on  factors  such 
as  chemical  composition,  concentration,  temperature,  electrode 
potential,  and  the  interaction  between  electrode  and  anions.  The 
literature  [68]  suggests  that  all  the  aromatic  compounds  interact 
via  their  aromatic  ring  parallel  to  the  Pt  surface,  the  benzene  and 
toluene  mainly  adsorbing  without  dissociation.  One  study  [63]  also 
indicated  that  the  adsorbed  toluene  can  be  electrochemically  oxi¬ 
dized  partially  to  higher  oxidation  state  intermediates  and  partially 
to  CO2. 

To  study  the  adsorption  behaviour  of  toluene  on  Pt  in  conditions 
similar  to  those  of  fuel  cell  operation,  a  toluene  adsorption  isotherm 
study  of  polycrystalline  platinum  electrode  was  conducted  [69]. 
It  was  found  that  the  Langmuir  type  of  adsorption  fits  well  the 
adsorption  isotherms  in  the  studied  temperature  ranges  (25-80  °C). 
The  estimated  free  energy  of  adsorption  revealed  toluene’s  strong 
affinity  towards  the  Pt  surface;  it  was  found  that  toluene  adsorbed 
mainly  as  a  cr-complex  (vertical  orientation)  at  low  concentration 
and  temperature,  while  the  TT-complex  predominated  at  high  con¬ 
centration  and  particularly  high  temperatures. 

Based  on  the  experimental  evidence,  we  propose  that  toluene 
adsorption  (reaction  (a-iv))  is  the  dominant  mechanism  for  air-side 
toluene  contamination  in  a  fuel  cell.  At  high  electrode  poten¬ 
tial,  toluene  electrochemical  oxidation  takes  place,  which  yields 
intermediates  and  carbon  dioxides  in  higher  oxidized  states  (as 
expressed  in  reactions  (a-vii)  and  (a-viii)). 


where  £Cgdl  is  the  porosity  of  the  cathode  gas  diffusion  layer  and 
^o2,air  is  the  bulk  oxygen  diffusivity  in  air. 

The  oxygen  concentration  in  the  cathode  catalyst  layer  is  calcu¬ 
lated  by  the  gas  law 


3.  Experimental  results 

The  testing  platform  design  for  fuel  cell  toluene  contamination 
has  been  described  in  another  paper  [70].  All  tests  were  con¬ 
ducted  on  a  Fideris  100  W  fuel  cell  test  station  with  a  single-cell 
hardware  purchased  from  Teledyne  (50  cm2  CH-50).  The  flow  field 
plates  were  designed  and  fabricated  in-house  using  single  ser¬ 
pentine  flow  channels  of  1.2  mm  width,  1.0  mm  channel  depth, 
and  1.0  mm  landing.  The  membrane  electrode  assembly  (MEA) 
has  an  active  area  of  50  cm2.  The  gas  diffusion  layer  (GDL),  pur¬ 
chased  from  SGL  Carbon  Group,  was  a  PTFE  (20wt%)  and  carbon 
black  impregnated  carbon  paper.  The  catalyst  coated  membrane 
(CCM)  (manufactured  by  Ion  Power)  was  made  of  Nation®  211 
membrane  with  0.4  mg  cm-2  Pt  loading  on  both  anode  and  cath¬ 
ode  sides.  A  fresh  MEA  was  employed  for  each  contamination 
test. 

In  all  fuel  cell  tests  both  with  and  without  toluene  present, 
the  voltage-current  polarization  curves  (steady-state  polarization 
curves)  were  collected  using  a  load  bank  controlled  in  a  constant- 
current  pattern.  The  relative  humidity  (RH)  of  the  fuel  cell  was 
controlled  at  80%  for  both  anode  and  cathode  sides,  and  the  cell 
temperature  and  backpressure  were  controlled  at  80  °C  and  30  psig, 
respectively. 

4.  Results  and  discussion 

4.1.  ORR  parameter  study 

In  order  to  obtain  ORR  parameters,  we  simulated  experimen¬ 
tal  baseline  data,  which  is  free  of  contaminants.  The  parameters 
employed  for  the  simulation  are  listed  in  Table  1,  and  the  simula¬ 
tion  and  experimental  polarization  curves  are  illustrated  in  Fig.  1 . 
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Table  1 

Parameter  values  used  in  polarization  fitting. 


F 

96485  C  mol”1 

R 

8.315  JK-1  mol-1 

T 

353  K 

»o2 

4 

Yc 

5000  cm2  cm-2 

rh 

OV 

n 

1 

n2 

1 

na3 

1 

Wap3 

0.5 

Wap3 

1 

0ip4 

0.5 

Q 

1 

OipS 

0.5 

z 

1 

f-CGDL 

3.53xl0-2  cm 

PCGC 

30psig 

A)2,air 

1.366  x  10-1  cm2  s”1 

£cgdl 

0.8 

Cm 

1.7273  x  10“3  mol  cm-3 

Table  2 

Estimated  parameter  values  based  on  the  baseline  polarization  curve. 

feu 

feib 

5.8  x  104  cm3  mol-1 

ka3b 

5.52  x  104  (cm3  )3  (mol-1  )3 

/<2f 

1.64  x  10-11  s_1  cm3  mol-1 

«2 

0.815 

Table  2  lists  the  estimations  of  these  parameters  from  the  simu¬ 
lation.  Our  estimated  rate  constant  for  the  rate-determining  step 
(reaction  (a-ii))  is  1.64  x  10-11.  This  number  is  quite  similar  to  the 
rate  constant  of  2.54  x  10-11,  which  is  derived  from  a  reported 
experimental  exchange  current  density  obtained  at  fuel  cell  operat- 


Fig.  1.  Comparison  of  calculated  and  experimental  polarization  curve  in  the  absence 
of  contaminant. 


ing  conditions  in  a  high  current  density  region,  with  a  temperature 
of  80  °C  [71].  The  estimated  symmetric  parameter  for  the  same 
reaction  (a2)  is  0.815. 

As  discussed  in  another  paper  [38],  feed  stream  contaminant 
concentration  at  the  catalyst  layer  is  a  function  of  both  inlet 
contaminant  concentration  and  current  density.  Furthermore,  the 
contaminant  adsorption  (desorption)  rate  constant  is  also  related 
to  electrode  potential.  We  can  treat  the  product  of  the  contaminant 
adsorption  (desorption)  rate  constant  and  the  contaminant  concen¬ 
tration  at  the  CCL  as  a  function  of  current  density  and  contaminant 
inlet  concentration  (/<Cp  ~/(Cp,/)),  where  Cp  is  the  contaminant 
concentration  in  the  cathode  channel.  With  parametric  study  of 
experimental  data  at  current  densities  of  0.2, 0.5, 0.75  and  1 A  cm-2 
and  contaminant  concentrations  of  1, 5,  and  10  ppm  (Fig.  2),  we  have 


100  150 

Time,  Min 


250 


Fig.  2.  Experimental  performance  data  and  modeling  results  at  different  toluene  concentrations  and  different  current  densities:  (a)  Jcen  =  0.2  A  cm-2 ;  (b)  Jceu  =  0.5  A  cm-2 ;  (c) 
at  /Ceii  =  0.75  A  cnrr2  and  (d)  /ceii  =  1.0  A  cm-2. 
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Fig.  3.  Comparison  of  simulated  and  experimental  steady-state  performance  under 
different  toluene  contaminations. 

obtained  a  relationship  between  contaminant  adsorption  (desorp¬ 
tion)  rate  constants,  toluene  concentration  at  the  CCL,  and  current 
density  that  enables  us  to  predict  cell  performance  and  surface  cov¬ 
erage.  Fig.  3  presents  a  comparison  of  steady-state  cell  performance 
under  different  toluene  concentration  levels. 

4.2.  Model  simulation 

With  the  developed  toluene  contamination  model,  we  simu¬ 
lated  transient  and  steady-state  cell  performance  (Figs.  4  and  5)  at 
different  current  densities  and  with  toluene  concentrations  at  ppb 
levels,  which  closely  resemble  normal  indoor  and  outdoor  toluene 
levels.  As  presented  in  Fig.  4,  at  a  current  density  of  0.2  A  cm-2 
the  voltage  drop  due  to  the  presence  of  750  ppb  toluene  can  be 
9  mV.  If  the  current  density  is  increased  to  0.5, 0.75  and  1.0  A  cm-2, 


Fig.  5.  Effects  of  toluene  contaminations  on  steady-state  fuel  cell  performance. 


the  same  amount  of  toluene  can  cause  corresponding  cell  perfor¬ 
mance  drops  up  to  16,  27,  and  48  mV,  respectively.  As  shown  in 
Fig.  4(d),  at  a  current  density  of  1.0  A  cm-2,  the  cell  voltage  drops 
due  to  toluene  concentrations  of  250,  500  and  750  ppb  are  37,  42 
and  48  mV,  respectively.  Fig.  5  shows  the  effects  of  toluene  contam¬ 
ination  on  steady-state  cell  performance  at  different  concentration 
levels.  Thus,  the  extent  to  which  toluene  contamination  affects  cell 
performance  depends  on  toluene  concentration  and  current  den¬ 
sity.  Based  on  this  model,  we  can  estimate  the  maximum  allowable 
toluene  concentration  in  order  to  limit  the  voltage  drop  to  a  certain 
range.  For  instance,  to  limit  the  contamination  potential  drop  to 
less  than  10  mV,  with  a  current  density  of  1.0  A  cm-2,  toluene  con¬ 
centration  should  be  less  than  100  ppb.  One  can  also  estimate  the 
degree  of  cell  performance  degradation  at  a  certain  contaminant 
level  and  current  density. 


Fig.  4.  Simulated  fuel  cell  transient  performance  at  different  current  densities  and  different  toluene  concentrations. 


Pt-02  surface  coverage  ( 0  )  ^  Pt-02  surface  coverage  ( 0 
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Fig.  8.  Oxygen-containing  intermediates  surface  coverages  at  different  current  densities  and  different  toluene  concentrations. 


Fig.  9.  Toluene  surface  coverages  at  different  current  densities  and  different  toluene  concentrations. 
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Fig.  10.  Surface  coverage  variation  with  overpotential  at  different  toluene  concentrations:  (a)  C°  =  Oppb;  (b)  C°  =  250  ppb;  (c)  C°  =  500  ppb;  (d)  C°  =  750  ppb. 


With  our  model,  we  investigated  the  effect  of  toluene  con¬ 
tamination  on  surface  coverage  at  different  current  densities.  The 
simulation  results  are  illustrated  in  Figs.  6-9.  Fig.  6  demonstrates 
the  effect  of  contaminant  level  and  current  density  on  Pt  surface 
coverage,  0pt.  Two  observations  can  be  drawn  from  this  figure:  (1) 
the  free  Pt  site  coverage  decreases  as  contamination  level  increases, 
indicating  that  toluene  contamination  blocks  the  free  Pt  sites;  (2) 
with  the  same  toluene  concentration,  increasing  current  density 
will  result  in  a  low  number  of  free  Pt  sites.  For  example,  at  a  current 
density  of  0.2  A  cm-2,  the  free  Pt  site  coverage  is  reduced  from  59% 
to  48%  as  the  toluene  concentration  increases  from  0  to  750  ppb, 
while  at  a  current  density  of  0.5  A  cm-2,  the  free  Pt  site  coverage 
changes  from  about  60%  to  about  40%.  The  same  observation  can 
be  made  of  Fig.  7,  which  simulates  the  surface  coverage  0Pt_o2  • 
Fig.  8  reveals  the  effects  of  current  density  and  contaminant  level  on 
oxygen-containing  intermediates  coverage,  0pt-o2H-  At  the  current 
densities  presented  in  this  figure,  the  oxygen-containing  interme¬ 
diates  coverage  is  small,  as  the  subsequent  reactions  are  much  faster 
than  its  formation.  Toluene  contamination  further  reduces  the  cov¬ 
erage,  and  the  extent  of  the  effect  is  greater  with  a  higher  toluene 
concentration.  With  respect  to  contaminant  surface  coverage,  Fig.  9 
illustrates  that  increasing  the  contamination  level  and  current  den¬ 
sity  will  in  both  cases  increase  the  contaminant  coverage,  0Pt-p. 
Depending  on  the  current  density,  with  750  ppb  toluene  the  0pt-p 
can  be  as  high  as  73%,  suggesting  a  significant  degradation  in  cell 
performance. 

We  also  examined  the  surface  coverage  of  all  the  species  at 
different  overpotentials  and  with  different  toluene  concentrations 
(Fig.  10).  When  contaminant  is  absent,  our  simulation  shows  that 
in  the  overpotential  range  from  0  to  0.23  V  the  catalyst  surface  is 
covered  mainly  (coverage  >50%)  by  oxygen-containing  intermedi¬ 
ates,  0pt— o2h-  From  our  model  simulation,  at  an  overpotential  of 
around  0.29  the  oxygen-containing  intermediate  coverage  is  about 
10%.  The  same  coverage  was  reported  as  well  at  a  cell  potential  of 


0.84  V  [40].  Furthermore,  at  an  overpotential  greater  than  0.26  V, 
the  free  Pt  sites  dominate  (i.e.,  coverage  >50%).  Oxygen  coverage 
(0pt-o2 )  is  also  generally  less  than  half  of  the  free  Pt  site  coverage. 

When  toluene  is  present,  we  observe  that  at  low  overpotentials 
(less  than  0.25  V),  the  toluene  coverage  (0pt-p)  is  small  regardless  of 
toluene  concentration.  This  is  consistent  with  experimental  obser¬ 
vation  [70]  that  toluene  contamination  has  an  insignificant  effect 
on  OCV.  Significant  toluene  coverage  (>10%)  is  mainly  observed  at 
large  cathode  overpotentials  (>0.32  V).  With  higher  toluene  levels, 
significant  toluene  coverage  will  appear  at  relatively  lower  overpo¬ 
tentials.  With  750  ppb  toluene,  10%  toluene  coverage  is  noticed  at 
a  cathode  overpotential  >0.30  V.  In  addition,  with  the  presence  of 
contaminant  we  see  a  narrowing  of  the  potential  window  in  which 
free  Pt  sites  dominate.  With  0  ppb  toluene,  this  free  Pt  site  window 
is  at  an  overpotential  >0.26  V;  with  250  ppb  toluene,  the  window 
is  >0.26  V  and  <0.38  V  (0.12  V  range);  with  500  ppb,  the  window 
is  >0.26  V  and  <0.37  V  (0.11  V  range);  and  with  750  ppb  toluene, 
the  window  is  >0.27  V  and  <0.36  V  (0.09  V  range).  Thus,  contam¬ 
ination  significantly  reduces  the  potential  range  in  which  free  Pt 
sites  dominate.  As  a  result,  the  fuel  cell’s  performance  is  adversely 
affected. 

4.3.  Model  limitation 

In  this  first  general  air-side  contamination  model,  we  have  con¬ 
sidered  oxygen  adsorption  on  a  single  Pt  reaction  active  site  to  avoid 
mathematical  challenges.  While  oxygen  adsorption  on  Pt  initially 
forms  an  end-on  configuration,  it  then  develops  a  bridged  configu¬ 
ration,  which  may  occupy  two  or  more  Pt  sites.  How  the  other  site 
model  fits  the  experimental  data  definitely  requires  more  investi¬ 
gation.  With  respect  to  a  toluene  contamination  model,  toluene 
adsorption  studies  indicate  that  the  adsorbate  configuration  is 
dependent  on  toluene  concentration,  temperature,  and  electrode 
potential.  Further  examination  of  the  number  of  adsorption  sites  is 
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underway  in  our  lab.  It  is  known  that  the  ORR  rate  depends  on  the 
Pt  crystalline  faces,  as  does  the  toluene  adsorption  rate.  This  rate 
difference  and  interplay  may  be  one  reason  why  overall,  our  simple 
adsorption  models  fit  the  experimental  data.  Furthermore,  in  real¬ 
ity  toluene  concentration  in  air  is  at  the  ppb  level,  so  experimental 
validation  at  the  ppb  level  is  in  progress  at  our  lab. 

5.  Summary 

We  successfully  developed  a  general  air-side  feed  stream  con¬ 
tamination  model,  which  has  been  validated  with  experimental 
data  on  fuel  cell  toluene  contamination  at  four  different  current 
densities  (0.2,  0.5,  0.75,  and  1 A  cm-2)  and  three  toluene  concen¬ 
trations  (1,  5,  and  10  ppm).  This  model  employs  the  associative 
oxygen  reduction  mechanism  and  considers  possible  contaminant 
reactions,  such  as  surface  adsorption,  competitive  adsorption,  and 
electrochemical  oxidation.  The  model  has  the  capability  to  simu¬ 
late  and  predict  both  transient  and  steady-state  cell  performance. 
Using  it  we  examined  toluene  contamination  effects  on  cell  per¬ 
formance.  Our  model  reveals  that  with  750  ppb  toluene,  a  48-mV 
performance  drop  at  1 A  cm-2  can  be  expected.  The  contamina¬ 
tion  level  and  current  density  determine  the  extent  of  performance 
degradation.  Our  model  simulation  indicates  that  without  contam¬ 
inant,  the  oxygen-containing  species  dominate  the  catalyst  surface 
at  overpotentials  <0.25  V.  At  overpotentials  >0.26  V,  free  Pt  sites 
dominate.  In  the  presence  of  toluene  contaminants,  the  surface 
coverage  is  not  affected  at  low  overpotentials  (<0.25  V).  Toluene 
contamination  dominates  at  larger  overpotentials  (>0.3  V).  With 
750  ppb  toluene,  20%  toluene  coverage  is  predicted  at  cathode  over¬ 
potentials  >35  V,  while  50%  coverage  is  predicted  at  overpotentials 
>0.40  V. 
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